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Abstract 

We present our full analysis of the two flavor Nambu-Jona-Lasinio model with 
SU{2) X SU{2) chiral symmetry on the four-dimensional hypercubic lattice with 
naive and Wilson fermions. We find that this model is an excellent toy field theory 
to investigate issues related to lattice QCD. We use the large N approximation to 
leading order in to obtain non perturbative analytical results over almost the 
whole parameter range. By using numerical simulations we estimate that the size 
of the corrections for most of the quantities we consider are small and in this 
way we strengthen the validity of the leading order large N calculations. We ob- 
tain results regarding the approach to the continuum chiral limit, the effects of the 
zero momentum fermionic modes on finite lattices and the scalar and pseudoscalar 
spectrum. 
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1 Introduction 



At present the low energy properties of QCD can only be studied numerically using lattice 
gauge theory. It is believed that in order to obtain reasonable results in reasonable time 
periods the computing power needed is of the order of one Teraflop. Although supercom- 
puters with such capacity may be built in the next few years it is still important to use 
other methods to get insights into the physics as well as into the behavior of the theory on 
the lattice. 

The Nambu-Jona-Lasinio (NJL) model was introduced before the discovery of quarks 
as a theory of nucleons interacting with a four-Fermi interaction. Today the fermionic fields 
of the model are reinterpreted as being those of the quarks. The most important feature 
of the model that qualifies it to describe some of the important low energy properties of 
QCD is that it possesses the same chiral symmetry as QCD and that this symmetry can 
be realized in the Goldstone mode. 

The NJL model can also be motivated by an argument found in The gauge field in 
the full theory of QCD develops a finite correlation length of the size of the inverse mass 
of the lightest glueball (~ 1550 MeV, see for example |@]). If we consider QCD on a lattice 
and integrate the high momentum fluctuations of the fermionic and gauge fields down to 
A < 1550 MeV then the lattice spacing will be of the order of the correlation length and 
we must then essentially have a theory of fermions with contact interactions and cutoff A. 
The resulting effective Lagrangian will maintain the original chiral symmetry but will of 
course be more complicated. If we further restrict our attention to energies much below the 
cutoff, naively speaking, it should be enough to keep in the Lagrangian the least irrelevant 
operator namely the four-Fermi dimension six operator. This is the NJL model. 

Unfortunately, by only keeping the four-Fermi operator, valuable information was lost 
and the model does not confine the quarks. Therefore, strictly speaking, it can not be a 
true effective field theory of QCD. Furthermore, if, for example, we want to study the a 
particle, which on phenomenological grounds is believed to have mass ~ 750 MeV, then 
the separation of scales is probably not large enough to justify the neglect of operators with 
dimension higher than six. Nevertheless, as mentioned above, the NJL model possesses the 
same chiral symmetry as QCD and it can realize this symmetry in the Goldstone mode. 
It is this feature that is most crucial in the understanding of the lightest hadrons and the 
reason for the successful quantitative predictions of the model. 

The NJL model has been studied extensively for various cases with continuum type 
regularizations. For a comprehensive review the reader is referred to [0] and references 
therein. Also for some recent work in the continuum see Furthermore the NJL model 
is a special case of Yukawa models that, under a different context, have been studied 
extensively with lattice regularization [Q. The model has also been studied on the lattice 
in connection with the possible equivalence of the top quark condensate with the Higgs 
field 0. In that work, however, the separation of scales is very large (the cutoff is of the 
order lO^'^ GeV), and it is therefore quite a different problem than the one considered here. 

In this paper we do not attempt to use the NJL model to make physical predictions. 
Our interest in the NJL model is purely qualitative and originates from the fact that 
its lattice version is an excellent toy model to investigate issues related to lattice QCD. 
This is the case primarily because it is an effective four dimensional fermionic theory that 
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has the same chiral symmetry as QCD, which symmetry in an appropriate phase may be 
reahzed in the Goldstone mode, and because one can use the large approximation to 
leading order in to obtain non perturbative analytical results over almost the whole 
parameter range. Also by using numerical simulations of the model one can estimate the 
size of the corrections and in this way strengthen the validity of the leading order 
large calculations. As it will be shown in this paper the 1/A^ corrections estimated in 
this way are small for most of the quantities we consider. This means that one has a good 
analytical handle on the model. Because of this a wealth of information can be extracted. 

For naive fermions, we calculate the scalar and pseudoscalar spectrum as a function of 
the cutoff and find agreement with general expectations. This is not perceived as a quan- 
titative prediction but rather as a first approximation to a physical picture. Furthermore, 
since the NJL model possesses the right symmetry of a possible embedding theory of the 
Higgs Sector it can therefore provide an interesting example of such a theory. We briefly 
discuss this possibility. 

One of the main purposes of our work involves the study of the Wilson formulation of 
lattice fermions. This formulation breaks the chiral symmetry of the model explicitly. The 
symmetry is expected to be restored as a result of the tuning of the Wilson parameter k 
to some critical value k^. Needless to say this is a very important issue for QCD and it 
would be very useful to clarify how it is realized in a simpler model. This involves, for 
example, issues such as what the continuum chiral limit on a lattice is and how does one 
define Kc. For Wilson fermions we also determine the effect of the heavy doublers on the 
sigma particle. We show that the heavy doublers raise the mass of the sigma to be of order 
cutoff. 

Another issue investigated in this paper involves the study of the effects of the zero 
modes on finite lattices. This is of particular importance to numerical simulations. We 
demonstrate that in certain cases these effects are large and can obscure the extrapolation 
to infinite volume physics. We also estimate the size of the 1/A^ corrections by comparing 
leading order large A^ results on finite size lattices with numerical results on the same size 
lattices and find that in most cases the 1/A^ corrections are small. This is of interest since, as 
mentioned above, the NJL model has been and is studied extensively as a phenomenological 
model using continuum type regularizations and among other methods the leading order 
large A^ approximation. 

In this paper we consider the two flavor (up and down) NJL model with SU (2) x SU (2) 
chiral symmetry and SU(N) color symmetry, with scalar and pseudoscalar couplings on 
the four dimensional hypercubic lattice. We consider both naive and Wilson fermions and 
we study the model using a large A^ expansion as well as a Hybrid Monte Carlo (HMC) 
numerical simulation. A shorter version of this work containing only the main results has 
appeared elsewhere 0. 

Throughout the paper we use a notation whereby small letters denote quantities in 
lattice units and capital letters denote quantities in physical (MeV or GeV) units. The 
lattice spacing is denoted by "a" . 

The paper is organized as follows. The model and its lattice version are described in 
Section 2. The large A^ analysis to leading order in large A^ is given in Section 3. Using 
asymptotic expansions we study in Section 3.1 the case of naive fermions on an infinite 
lattice and present results regarding the scalar-pseudoscalar sector of QCD and briefly 
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discuss their possible relevance to the Higgs sector. In Section 3.2 we study the case 
of Wilson fermions on an infinite lattice using asymptotic expansions and present results 
regarding the approach to the chiral continuum limit as well as the effect of the doublers on 
the scalar particle self energy. The numerical and large work on finite lattices is described 
in Section 4. We present there, in Section 4.1, a comparison of our numerical results with 
the large results obtained on same size lattices and obtain an estimate of the size of the 
1/N corrections. In Section 4.2 we examine the effects of the zero momentum quark modes 
on finite lattices in connection to the inversion time of the Conjugate Gradient algorithm 
used in the HMC, and more importantly in connection to finite size effects. Furthermore 
after the preprint version of this work appeared in the computer list "hep-lat" and during 
the refereeing process of this paper another group [|l^ presented complementary results. 
In particular, for the case of Wilson fermions, the authors showed, using the leading order 
large N approximation, that there is a phase where the remaining flavor-parity symmetry 
breaks spontaneously generating a non zero vacuum expectation value for the pion fleld. 
Their work was done at inflnite volume using numerical integration and among other things 
it was shown that the phase line is a line at which all three pions are massless. This line 
was plotted for all effective quark masses. In this work we always stay on the symmetric 
side of this line where the vacuum expectation value of the pion fleld is zero. At inflnite 
volume, using asymptotic expansions, we had calculated the part of this line (flg. 2a) 
that corresponds to small effective quark masses since this is the region where continuum 
physics is extracted. Our calculations, up to errors relevant to the approximations used, 
are in agreement. However, on a flnite volume the comparison is not as straight forward 
and it deserves special attention. For this reason section 4.3 was added where the zero pion 
mass line on a flnite volume is discussed in detail. Finally, a short summary and conclusion 
is given in Section 5. 



2 The model 



The Lagrangian density in Minkowski space and in continuum notation is: 

C = ^(z^ - mo)^ + ^ + (^i75T^; 



(1) 



In the above expression all indices have been suppressed. The fermionic fleld \1/ is a flavor 
SU (2) doublet and a color SU{N) A^-column vector. The Lagrangian is diagonal in color, 
in contrast with the full QCD Lagrangian which is diagonal in flavor, r = {ti,T2,T3} are 
the three isospin Pauli matrices, ^ = 7'^^^, and mo is the bare quark mass (if mo 7^ 
the chiral symmetry is explicitly broken). To obtain a Lagrangian that is quadratic in 
the fermionic flelds we introduce the scalar auxiliary fleld a and the three pseudoscalar 
auxiliary flelds tt = {tti, 7r2, tts}. Using the functional identity M 



exp \ I d'^x {^^f + (^i75T^^2 



j [dad-n] exp |i j d^x -^(cr + i'jr^r • tt)^ - n//3i(a^ + tt^) | 



(2) 
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the Lagrangian density becomes 

C = M^-n//3i(a2 + 7r2) 

M = — rriQ — a — i'y^T • tt . (3) 

Here n/ = 2 is the number of flavors and Pi = ^^^^^ . The fermionic fields can now be 
integrated and the resulting partition function is 

Z = J [dadn] [det Mfe-'^'f^' / d^x{a^+7z^) _ 

Notice that there is no explicit kinetic energy term for the a and tt fields. As we will show 
soon this term is part of det M. 

In the general case the fermionic determinant has a phase that is related to the the 



chiral anomaly and the Wess-Zumino term [11]. In our work because the flavor space is 



restricted to SU (2) the phase is not present and we will write our action in a symmetric 
fashion with regards to M and . Going to Euclidean space and appropriately discretizing 
the above Lagrangian we obtain the model on the Euclidean Hypercubic lattice. On the 
lattice, as it is well known, we have species doubling. The doubling in the NJL model will 
be interpreted as a doubling of the color degrees of freedom. To treat this problem we 
add to the Lagrangian density an irrelevant operator (Wilson term) of the form ^^t^^^; 
"a" being the lattice spacing and "r" a constant. We consider the r = case where no 
effort is made to remove the doublers (naive fermions) and also the r ^ case where 
the doubler masses are raised to the cutoff (Wilson fermions) and the chiral symmetry is 
explicitly broken. With these considerations and after appropriate scaling of the flelds and 
couplings, so that only dimensionless quantities appear, we obtain: 

[d^d'^dad7T]e~^ 

{N/2 
1=1 ^ ■' 

with 7^ hermitian. This partition function can be cast into two different forms. The flrst 
can be studied using a large N expansion and is given by: 



[dad7r]e' 



Si = iv|-lTr(logM)-lTr(logMt) + n;/3i^(a^ + 7r^)| (6) 

with Pi = Pi/N and the trace taken over space, spin, flavor, and color. The second is 
appropriate for numerical simulations and is given by: 



j [dXdXUad7z]e-^^ 

'N/2 

' it, 



= Y.\Y.{^x{M^M)l^Xl]^n^Pi{al^^\)b^y\ (7) 
with X being pseudofermionic flelds. 
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3 Large 



We perform a standard large expansion with the action of eq. ^. The large ap- 
proximation will be reasonable for as long as (3i is of order one. We assume a translation 
invariant saddle and small fluctuations around it. For the case of naive fermions using the 
chiral symmetry we rotate the fields so that the saddle field configuration lies along the a 
direction. For the case of Wilson fermions the chiral symmetry is broken and the saddle 
field configuration lies along the a direction. As was shown in |T^, the remaining parity- 



flavor symmetry can spontaneously break generating a non zero vacuum expectation value 
for the pion field. However, in this work we will always stay in the parity-flavor symmetric 
phase where the pion field has zero expectation value. We have: 

a(x)=a, + — U, 7r(x) = — U (8) 



In momentum space the inverse quark propagator at the saddle is: 

^s(p) = ^Xl^M^i^^'M + ^(4 - XI'^osPm) + "^0 + • (9) 

We identify 

mo + cr^ = niq (10) 

as the quark mass. This definition is valid close to the continuum limit where is small. 
It must be pointed out that is the constituent quark mass since all high energy gluonic 
degrees of freedom of QCD have been integrated out. In physical units the quark mass will 
be taken to be equal to one third the proton mass. 

We expand around the saddle (expansion in ^^^)- The zero order gives the 

effective potential 

^ = -2ni ^\og\g{j>.m,)\^nffi^al (11) 

where = /^^^ for an infinite lattice (B denotes the hypercubic lattice Brillouin 
zone) and = zfTI ^{ni,n2,n3,n4} with ni,n2,n3 G [0, - 1], G [0, - 1] and p = 
/2mii 2^ 2mi3 2mi4-) fgj, ^ finite lattice with spatial extend and temporal extent Lt- In 

L/x J^x J^t 

the above equation g{p,niq) is defined as 

g{p,mq) = ^ sin^ py + [rw{p) + mqf = p'^ + m'^ + 0{p'^mg) 

V 

2 

w{p) = A-Y^cosp, = ^ + Oip') . (12) 

The saddle point equations are obtained at the point where the linear term vanishes. 

/3i f as + mo + rw{p) 

- / — J — ^ = 13 

2 Jp g{p,mo + as) 
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The (J and tt propagators Ga{q) and G.,r{q) can be obtained from the second order term 
of the expansion around the saddle | J^6a{q)G~^6a{—q) + \ j^5Tz{q)G~^5iT{—q) with: 



G-/M) = 4ni 



'/3 



h/a{q) 



Y.uSi\i{K + ^)si\i{K-^) +/- rw{k + 



rrir. 



rw{k- 



(14) 



(15) 



g{k + ^,mq)g{k-^,mq) 
We define the pion wave function renormahzation constant and pion mass m,r from: 

hmG^'iq) = Z;\q' + ml) . (16) 



With this definition we find using eq. 14 and 15: 



ry-l 2 

Z„ mt 



Uf f Em cos^ kf, + sin^ k^. 



An, 



7?1 



g{k,mqy 



k g{k,mq) 



(17) 
(18) 



Notice that is positive and therefore does not create local stability problems (the prop- 
agator is positive). Also notice that with this definition of the pion mass, will be the 
true pole of G^r only if m^r is small. For the cases we are interested in this will always be 
true. Unfortunately we can not define the a mass mg- and width 70- in a similar way since 
TJirj will not be small in general. The proper definition in this case is: 



G-\q = {im„ + ^MQ]) 







(19) 



and we can not obtain closed form expressions for ma and 70-. 

On a finite lattice the finite momentum sums can be calculated with the aid of a 
computer and the above quantities except 70- can be exactly determined. This will be done 
in Section 4.1 when we will compare the numerical results obtained on a finite lattice with 
the large N predictions for the same lattice size. For an infinite lattice the momentum 
integrals can be evaluated numerically but this will not serve our purposes since we are 
interested in obtaining analytical expressions. For this reason we resort into performing 
asymptotic expansions in small m^, m^r, mo-, and g, considering the logarithms as being of 
order zero. The asymptotic expansions will be given separately for the r = and r 7^ 
cases in the next two sections. 



3.1 Naive fermions at infinite volume, r = 

In this case we have species doubling. The doubling is interpreted as a "doubling" of the 



color degrees of freedom. That this interpretation is appropriate can be seen from eq. ITT 
For r = the function g{p^ m^) of eq. is periodic with period vr and not 2tt. This means 
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that the integral in eq. 11 sphts into 16 equal pieces and therefore the effective potential 
will be made up from 16A^ copies. Therefore we set the number of colors N^. to: 



iV, = 16A^ . 

Some of the typical integrals we will encounter are: 

J„(mg) 



(20) 



k g{k,mq 



(21) 



Ik g{k + ^,mg)g{k-^,mg) 

Simple trigonometric relations relate these integrals with the corresponding integrals that 
arise for Bose particles. We get 



JJm) 



16{4"-2j^(2m)} 



/o(g,m) = 16/o^(2g,2m) 



(22) 



where and Iq are defined as in eq. ^ but with g{k, m) = 2 [1 — cos{ky)] + rn? the 
inverse propagator of a Bosonic particle. The asymptotic expansions of and have 



been evaluated to leading order in nig with very accurately determined coefficients in [|12 
Appendix B. The leading order term of Iq is universal except for a lattice constant. For 



example, it has been evaluated for the F4 lattice in [T^. The lattice constant needed also 
appears in the leading order term of the asymptotic expansion of and for the hypercubic 



lattice can be taken from [Q. For the convenience of the reader we present these results 
below: 



Jf(m) 



Tq + rn?{ri + Si logm^) + 0(4) 
1 d 



n dm 



2 " 



Ioiq,m) = -si f rfslog[m2 + s(l-s)g^]-(ri + Si) + 0(2) 
Jo 



(23) 



with 



ri 

Sl 



0.154 933 390 

-0.030 345 755 
1 



167r2 



(24) 



Because r = 0, if we set mo = the chiral symmetry is not explicitly broken. In that 
case there are two solutions to the saddle point equation 0, namely o"s = and cxs 7^ 0. For 
(3i > f3i^, the dominating saddle is cr^, = and we are in a chirally symmetric phase with 
massive a and tt fields and massless quarks nig = = 0. For /3i < j3i^ the dominating 
saddle is cr^ 7^ and we are in a phase with spontaneously broken chiral symmetry. The 
pions are the Goldstone bosons and become massless. The sigma is the massive mode and 
the quarks acquire a dynamically generated mass rrig = cr^ 7^ 0. The critical value of jSi is: 



= 2Ji(0) = 8ro 



(25) 
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All of our analysis will be done in the broken phase. 
The asymptotic expansion of G~^{q) is: 



G-\q) = Z-^ml + Anjq^ |l/o(g, m,) - ^Ji(O)} + 0(4) 



(26) 



where only the leading order term of Z~^rn^ and Jq is to be kept. The coefficient of the 
term will be equal to Z^, for g — > 0. The pion mass mass should satisfy G~^{q = 
{im,r, 0, 0, 0}) = 0. To this order in the asymptotic expansion this will be true only if the 
coefficient of the term for q = {zm^, 0,0,0} is very close to Zt,. From the asymptotic 
expansion of Io{q, rriq) (eq. |^) we see that this will happen for m^s(l — s) ^ where 
s G [0, 1]. Therefore for this definition of to be valid we must not only demand that 
m-Tr is small, as mentioned in the previous section, but also that: 



(27) 



For all the cases we will be interested in this will be satisfied. For example with M„ 



310 MeV and = 140 MeV we get ^ = 0.054. 



Equation 17 can be rewritten in the form: 



Z. 



-1 _ % 



(4 + ml)J2{mq) - Ji(mg) 



(2^ 



Using eq. 22, 23 and 28 we find the leading order term of the asymptotic expansion of Z^: 



Z^^ = 16[2;o - 2;ilogm^ 



zo = -2nf(— + ri + si + silog4), for rif = 2 = 0.022 204 130 
16 

Zi = 2nfSi . 



(29) 



As mentioned in the previous section Z^^ is always positive and does not create local stability 
problems. From the above equation we find that will become negative if rriq > exp ^ ~ 
1.55. This provides us with a point above which the small approximation is certainly 
not valid. In our analysis we will never need to be close to this point. However, the above 
equation has a much more important consequence. It can be shown ^ that the pion decay 
constant fn is related to rriq and Z^, via: 



Nm^Z; 



(30) 



Using eq. BG, 2£ and 20 we get: 



J 'I 



= Nc[zo - zilogm 



(31) 



Then with Nr 



Mqtt, 



A 



- we 



310 MeV ""-"^ ""1 

find A = 1150 MeV. This is consistent with the expectation that the cutoff of the theory 
should be ^ 1550 MeV (the mass of the lightest glueball 0). 
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Although, as mentioned in the introduction, it is not the aim of this work to study the 
possibihty of the Higgs being a top quark condensate, we point out that similar consider- 
ations as in the previous paragraph can be used to estimate the cutoff of such a theory. 
With = 3,nf = 1, = 246 GeV, Mg = 170 GeV we obtain A = 3 x 10^^ GeV. This 
is consistent with the much more complete and detailed analysis of M\. 



Equation 18 can be rewritten in the form: 



ml = AnjZ^ — — — Ji(m,) (32) 
mg — mo 



where the saddle point equation [T^ was used to eliminate Pi. Using eq. ^ and we 
find the leading order term of the asymptotic expansion of mj^: 

{mg - mo) [zo - zi logm^J 

Notice that when the chiral symmetry is not explicitly broken, namely mo = 0, then in the 
broken phase where mg — mo = ag ^ 0, m^^ = as it should. 



Next we concentrate on m^^ and 7o-. Using the saddle point equation 13 to eliminate 



(3i and eq. |2^ the leading order term of the asymptotic expansion of the inverse sigma 



propagator of eq. |TJ is: 



16 



^N,[zo - z, logm^lX + 4r^/ {^(g' + 4m^^)/o^(2g, 2m,) - ^g'} + 0(4) (34) 



where only the leading order term of /(f given in eq. |23| is to be kept. The leading order 
term of /(f contains a log(— 1) that we take to be equal to +i7r (taking it to be equal to 
— iTT results to a negative width). Using the above equation and the definition ^ we can 
calculate m^. and 7o-. Because we can not get a closed form expression we solve numerically 
and the result is shown in Figure 1. 

From Figure 1 we find that if we set the quark mass to one third the proton mass 
Mg = 310 MeV then = 726 MeV, = 135 MeV, and A = vr/a = 1150 MeV. 
is consistent with phenomenological expectations and as mentioned earlier A is consistent 
with the expectation that the cutoff should be close and below the mass of the lightest 
glueball (1550 MeV). The width however is underestimated. The reason is traced to the 
fact that to leading order in large N the width receives contributions only from the quark 
bubble and not from the pion bubble because the pion bubble is of order 1/A^. Because 
the phase space available for the a to decay to two quarks is much smaller than the phase 
space to decay to two pions, the pion loop contribution, although of order 1/A^, is probably 
more important than the quark loop contribution. 

The above result can also be used to make an interesting observation. If the Higgs sector 
is the low energy effective field theory of a NJL model with exactly the same parameters 
as the low energy QCD except for M^^ = and = 246 GeV, then we find M„ = 
1915 GeV. This corresponds to = 2 where one would expect very large deviations from 
the low energy behavior of scattering cross sections. Although we have not calculated these 
deviations the value of the width serves as an indication of their size. In a way, departure 
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from low energy behavior will be signaled by an increasing width of the a to two quark 
decay. At mg. 2 the width is already fairly large. 












m 



O" 



Figure 1. r = 0, A^^ = 3, = 93 MeV, and = 140 MeV. From top to bottom 
the lines correspond to the a mass, quark mass and a width calculated in the large N 
approximation to leading order in rrig. The vertical line denotes the point where the quark 
mass is equal to 310 MeV. 

As a final comment notice that close to the continuum limit and by setting the renor- 
malization point of the integral /(f in eq. at g = we obtain 



7a 



ml + 







(35) 



where we have used the definition of eq. |T9|. This is the result obtained in the literature 
with continuum type regularizations (see for example |^). This result, because of the off- 
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shell renormalization point, is only approximate and also neglects the contribution of the 
quarks to 70-. 



3.2 Wilson fermions at infinite volume, r 7^ 

In this case the doublers have been removed by raising their masses to the cutoff. The 
chiral symmetry has been explicitly broken by the Wilson term and the pions are massive. 
The number of colors is: 

N, = N . (36) 
Some of the typical integrals we will encounter are: 



k g{k,mq 



Io{q,mg) = I ^- — . (37) 

Jk g{k + ^, mq)g[k-^, nig) 

Unfortunately there are no simple trigonometric relations that relate these integrals with 
the corresponding integrals that arise for Bose particles as in the r = case. We will need 
the asymptotic expansions of Jo,i, Ji,i to order and the leading order term of Jo,2 • We 
have calculated them and the result is given in Appendix A. We will also need the leading 
order asymptotic expansion of the Iq integral. This, up to the lattice constant ri that needs 
to be calculated and is given in Appendix A, is the same as the leading order term of J^f 
in eq. |2^. Some of the lattice constants that we will need can be parametrized as: 

fln.m = Jn,m{0) m - U < 1 . (38) 

In our analysis only few of them appear and we have calculated them for r = 1 and r = 0.1. 
They are given in Appendix A together with the constants 

= lib ■ '^«) 

where 6 is the step function. 

The asymptotic expansion of G~^{q) is: 

G-\q) = Z-'ml + Anjq' |i/o(g, m,) + ^^^^«o,i - ^^^^^^^2,2} + 0(3) (40) 

where only the leading order term of Z'^irt^ and Jq is to be kept. The same restrictions 
on as in eq. ^ are needed. 

The leading order term of the asymptotic expansion of is: 



^ = [zq- zi log 



Zq = —2nf 



1 — r^(l — r^) 
ri + si^ — ao,i 02,2 



for n/ = 2 ,r = 1, zq = 0.0223 



zi = 2nfSi . (41) 
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Notice that for r = 1, zq has a value that is very close to the one of the r = case. Because 



of that and since zi is universal all the discussion relating to eq. |29 
here. 



and 31 is also valid 



Equation 18 can be rewritten in the form: 



mo 



[fnoJo,i{mg) +rJi,i(mg)] 



(42) 



where the saddle point equation was used to eliminate The leading order term of 
the asymptotic expansion of rrij^ is: 



Ant 



[rriq - mo) [zq - Zi logm^] 



+ 0(2) . 



(43) 



We are now in a position to investigate the approach to the continuum chiral limit. The 
theory has two adjustable bare parameters, (3i and mo. The parameter mo is related to the 
more familiar hoping parameter k (often used in QCD) by: 



K 



8r + 2mo 



(44) 



The two bare parameters control rriq and m^ through equations |1^ and The following 
statements can be made: 

A The rriq = line where the continuum limit is retrieved is defined by: 



mo 

B For any rriq, mo can be adjusted so that m^^ = 0. The m^ = line is given by: 

moJo,i{mq) + r Ji_i(mg) = 



(45) 



(46) 



As was shown in [|T2] this line separates the flavor-parity symmetric phase from the flavor- 
parity broken phase. In this work we always stay in the symmetric phase. 
C The continuum chiral limit is obtained when both rriq and m^ go to zero. This point is: 



^0 ^ , = —'T , Pi ^ , = 

00,1 ^0 



(47) 



and all physics should be extracted at the vicinity of this point. 

Using the asymptotic expansions of eq. |1^ and ^ to order m^, we plot for small rriq 
and m,r the constant rriq and constant m^ lines for r = 1 in Figure 2a and for r = 0.1 in 
Figure 2b. The two Figures are qualitatively the same with the rrtq = 0, m^ = point 
shifted toward larger k and Pi. In fact as r — > it can be easily shown that j3i — > N/3i^, 
where f5i^ is the critical value of the r = case (see eq. |2S|), and k — > oo. From eq. ^ 
we find that as ruq also m^^ ~ 1/logm^ — 0. This implies that the rUq = line is 
also an m^^ = line. This behavior is apparent in Figure 2a and 2b. However, this is an 
artifact of our approximation in defining mT^. As mentioned earlier (see the comment after 
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eq. 1^ ) our definition of will be valid only if -C which can not be satisfied on 
the rrig = line except on the one point where it intersects the = line. The constant 

lines are therefore valid only in the regions where -C 4m^. To get the full constant 

lines a more detailed analysis would be necessary 



0.45 



0.40 



0.35 



0.30 

0.25 0.30 0.35 0.40 0.45 

1 

Figure 2a. The (3i, k = 1/(8 + 2mo) plane for r = 1.0, = 2. The solid lines are constant 
171^^ lines. From top to bottom they correspond to m^r = 0,0.1,0.2,0.3. The dotted lines 
are constant rrig lines. From right to left they correspond to rrig = —0.2, —0.1, 0.0, 0.1, 0.2, 
0.3. The Trig = 0,m^ = point is located at Achira! ~ 0.3416, Hchirai = 0.3994. 

From these Figures we see that if for a fixed f3i we were to change k from smaller to 
larger values (as is often done in QCD with dynamical Wilson fermions) then if j3i < /^i^^^j^^i 
we would reach the m,r = limit before we reach the continuum limit. On the other hand 
if /3i > /5i ^ , we would reach the continuum limit before we reach the = limit. As 

' ' ^chzral " 

mentioned in "C" above, there is only one point in the k plane where we can obtain a 
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continuum chiral limit. 



Using the saddle point equation |13] the chiral condensate is found to be: 



-2/?i(mg - mo) = -2(3ias 



(4J 




and it is not zero when irtq = and/or m.,^ = 0. Therefore it can not serve as an order 
parameter for Wilson fermions. On the other hand when r = and rriQ = (\I'\E') is an 
order parameter because cr^ is an order parameter. 

Next we concentrate on rricr and 7o-. The leading order term of the asymptotic expansion 
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of the inverse sigma propagator of eq. [T3 is: 



G-\q) = -^iV,[;^o-^ilogmJ]X 



qi q 



1 fr^ - 1) 

+ ^nf{-{q^ + Aml)Ia{q,m^) + q^ ^ ^ 



ao,i + -Ro + mqRi + mqi?2 + g^-Rs 



+ 0(3) 



(49) 



where only the leading order term of Iq is to be kept. Rq, Ri, R2, and R3 are lattice 
constants defined below. 



Ro 
Ri 
R2 



4rai,2 - 8r^a3,3 
24r^a4.4 - 20r2a2 : 



r2-l 



16 



4 Jkg^{k,0) 



\ sin^ Igik, 0) - 2rV{k) - 2w{k) cos(A;^ 



(50) 



i?o, -Ri, R2 can be calculated from the a„_m's. -R3 has to be calculated separately and is 
given in Appendix A. Notice that Rq, Ri, R2 and -R3 go to zero for vanishing r. 

It is immediately apparent that the Rq and nigRi terms do not scale appropriately. As 
a result 

Ro (51) 



and therefore Mo- is of order cutoff. For the NJL model at large we can trace the 
reason for this phenomenon and offer an exact answer. Such a phenomenon may also be 
responsible for the difficulty in observing a a particle in numerical simulations of QCD with 
dynamical Wilson fermions. 

First we realize that although the Wilson term has raised the doubler masses to the 
cutoff, the doublers have not disappeared and they can possibly contribute through vacuum 
polarization effects. We must then try to separate their contribution. Towards this end 



consider the defining equation of the a propagator eq. |T4|. The /?i/2 term contributes to 



the term of eq. 



F2 



j. This term scales "correctly" and it will not concern us. We 

focus on Jo of eq. |15[ I^{q) ~ Re Tr J^M-^{k-^)M-^{k + l) where M^p) is the quark 
propagator at the saddle point and is given in eq. ^. This integral corresponds to a bubble 
integral with two external legs a{q) and cr(— g) and two quarks flowing in the bubble with 
momenta k — q/2 and k + q/2. Now, each species has momentum that belongs in a section 
Bi, z = 1, 2, . . . 16 of the Brillouin zone B with extent 7i/2 and — 7r/2 from the origin in 
each direction (see the first two columns of Table 1). For small q both quarks will have 
momentum around k and we can then separate the contributions of each species by splitting 
the integral into the 16 regions B^. This way we can isolate the contribution to each term in 
/o from the propagating quark and the 15 doublers. We will only discuss the contributions 
to the Ro term since ~ Ro, but similar arguments hold for the -Ri, R2 and R^ terms. 
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From the definitions in 38, 3? and splitting the integral into 16 regions as described 
in the previous paragraph we obtain: 



where we have explicitly denoted the r dependence of g{k,0). It is immediately apparent 
that the contribution, Rq^, form the Bi region where the propagating quark lives is not 
zero. Therefore we can not attribute the whole Rq term to the doublers. 

The contribution of the propagating quark comes from the high momentum section of 
the Bi region. To see this consider a spherical section Ai ~ centered around the origin 
(0,0,0,0) of El. We have: 

Since W{k) ~ k'^ and gr{k, 0) ~ k'^ (see eq. |T2[) we find g^(ko)2 ~ therefore 

The Wilson term has not only raised the masses of the doublers but has also changed 
the high frequency behavior of the propagating quark. In a bubble integral this change is 
visible. 

We could contrive a Wilson term that will raise the doubler masses but not change 
the high frequency behavior of the propagating quark. This can be done by introducing a 
momentum dependent r such that: 

r{p) = a p E Bi 

r{p) = r iipeB2, B^, . . . Bie . (55) 



Then 



Rn 



Ib 9r{k,0y f^^JB, gr{k,Oy 

and the a mass mg. ~ Rq will be composed entirely from contributions due to the doublers. 

The above discussion can become more transparent for r <^ 1. The masses of the 16 
species are the roots of sin^ k^, + r^(4 — cos A^j^)^. For small r they are given in Table 
1. 

We calculate the leading order term of the asymptotic expansion of Rq in rrii and find: 
Ro = 2r' f "^^L -^flm^ h + ^1 + ^1 logm^j + 0(4) (57) 

JB, Qr-=n{k. r ~f L J 



1=1 



with ri and Si given in eq. The first term comes because of the change of the high 
frequency behavior of the propagator of each species. The second term clearly shows that 
can be written as a sum over the doubler masses squared. 
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4 Numerical and large work on finite lattices 

We simulate the action in eq. |^ for N = 2 using a standard Hybrid Monte Carlo (HMC) 



algorithm [Q. We use the conjugate gradient (CG) algorithm without preconditioning to 
invert the matrix M'^M and the leap-frog algorithm to integrate the equations of motion. 
We measure expectation values of operators that involve the a and tt fields as well as 
operators that involve the pseudofermionic fields or traces of appropriate combinations of 
the matrix M. 

The reader will realize that the matrix M of eq. ^ of the NJL model does not provide 
any significant advantages as far as inversion time is concerned over the matrix M of full 
QCD. The CG will be as time consuming as in full QCD. The advantage comes because of 
another reason. Since the cutoff of the NJL model is at < 1550 MeV and since we will want 
to look at energies around the pion mass we are dealing with a ratio of scales ~ A/M^^ ^ 10. 
Therefore we should expect that lattices of size 16^ should be large enough for this purpose. 
These size lattices can be simulated in 1/4 of the 64K processor CM-2 supercomputer at 
SCRI in a reasonable amount of time. Although this was the original justification for 
performing the numerical simulation it turned out that the numerical results had a more 
important consequence. As we will demonstrate in Section 4.1 the numerical results on a 
given finite size lattice are in good agreement with the leading order large N results on the 
same size lattice, indicating that the 1/A^ corrections are small for the quantities we were 
able to measure. 

The large N results on finite size lattices are obtained by explicitly performing the 
four-dimensional finite momentum sums on a workstation. 
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Some of the typical parameters of our numerical simulation are: 

* Trajectory length r = 1 

* Step size dr ~ 0.02 - 0.05 

* CG residue 10"^ 

* CG iterations 20-160 

* Acceptance rate > 90% 

* Mesurements per "point" ~ 100 

* Autocorelations 2 — 10 

* Time for one CG iteration for a 16^ lattice on 16 K processors of the SCRI CM-2 
~ 1.3 sec 




Q_Q, I \ \ \ \ \ \ \ \ \ \ ^ \ \ \ \ \ \ 1 

1.5 2.0 2.5 3.0 

^1 

Figure 3. r = 0, mo = 0, = 2, = 8, = 16. The diamonds are the values of rrig 
and the crosses are the values of m'^ = {—N < > /{2f3i)} + ttlq from the numerical 
simulation. Because of the functional identity, eq. |, the two quantities are expected to be 
equal. The solid line is the large A^ prediction on the same size lattice. The dotted line 
denotes the infinite volume f3i^ from the large A^ calculation. For f3i < f3i^ the model is in 
the broken phase. 
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12 3 
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Figure 4. Same as in Figure 3 but for the U{2) x U{2) NJL and for r = 1, = 8. Notice 
that the diamonds and crosses are almost identical. 

4.1 Numerical results and large N 

In this section we present our numerical results for N = 2 and we compare them with the 
large N predictions (scaled to = 2 when necessary) on same size lattices. 

To check for consistency of the fermionic and auxiliary fields we plot rriq and m'^ = 
{—N < ^'\E' > /(2/5i)} + mo vs. jSi as determined from the numerical simulation. Because 
of the functional identity, eq. ^ we should have rriq = m'^. This relation is satisfied nicely. 

The agreement with large of dynamically determined quantities vs. bare quantities 
is quite good and helps us to get oriented in the bare parameter space. This can be seen 
for various values of the parameters in Figures 3, 4, 5, and 7a. 
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The important comparison with large that will help us get a feel for the size of the 
corrections comes from comparisons of dynamically determined quantities vs. other 
dynamically determined quantities. In particular we exchange one of the bare parameters 
for {a). 
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ITIq 

Figure 5. r = l, /3i = 2.5,A^ = 2. The diamonds are the values of rriq from the numerical 
simulation for L^. = 8, = 16. The solid line is the large N result on the same size lattice. 

In Figure 6 we present the a and tt propagators in momentum space for ten small 
momenta and r = 0. It is from this figure that we would have to extract Z^^. As it can 
be seen the large prediction for the same (a) as the one measured in the simulation is 
in good agreement with the numerical results. The large predictions in Figures 6a, 6b, 
6c and 6d "fit" the numerical results with per degree of freedom 0.38, 0.67, 0.32, 0.42 
respectively. This means that the determination of = ^ ^ function of (cr) has 

small 1/A^ corrections. 
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Figure 6. a) The a propagator in momentum space for 10 small momenta with r = 0, 



mo = 0, pi^jc = 2.2, (or) = as = 0.4840, N = 2, = 16, Lt = 16. The crosses are 
the values from the numerical simulation. The diamonds are the large N results and "fit" 
the numerical results with a x^/d.o.f. = 0.38. b) Same as in 6a but for Pij^^f^ = 2.4, 
(a) = as = 0.35. The large results "fit" the numerical results with a /d.o.f. = 0.67. 
c) Same as in Figure 6a but for the pion propagator. The large N results "fit" the numerical 
results with a x^/d.o.f. = 0.32. d) Same as in 6b but for the pion propagator. The large 
N results "fit" the numerical results with a x^/d-o.f. = 0.42. 

Another dynamically determined quantity that agrees well with the large prediction 
when plotted vs. (a) is the pion mass m7r. This plot is shown in Figure 7b. This figure 
suggests that the 1/A^ corrections to rriT^ are fairly small. In this figure, as well as in figure 
7a, was not calculated from the definition |16| but as the imaginary pole of G-,, (see eq. 
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0.125 0.130 0.135 0.140 

/c 

Figure 7a. pi = 2.5, r ^ 1, N ^ 2, = 8, Lt ^ 16. The crosses are the MC data. The 
sohd hne is the large N prediction on the same size lattice. 

It should be noted that the good agreement of the large results with the numerical 
simulations is not only present for naive fermions where the number of species is 16 times 
larger and therefore one would have naively expected the leading order large N expansion 
to be a good approximation. It is also present for Wilson fermions (r 7^ 0) as it can be 
seen from Figures 4, 5, 7a and 7b. 

Finally, as discussed in detail in Section 3.2, for r 7^ the a mass is of order cutoff 
and therefore very heavy to be able to measure from the decay of the a — a correlation 
function. However, for r = one would expect to be able to measure rria- As we will 
discuss in the next section this is not possible with the lattice sizes accessible to us. This 
is unfortunate since m„ is another very important quantity. However, we expect the size of 
the 1/N corrections of mo- to be similar to the ones of and therefore fairly small. Also, 
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measurements of the sigma width were not performed, but, as discussed in Section 3.1, we 
expect the corrections to the width to be fairly large. 
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Figure 7b. Same as in Figure 7a but as a function of (a). 



4.2 The zero momentum mode of the quarks 

To leading order at large the matrix M'^M of eq. |^ is diagonal in momentum, spin, 
flavor, and color spaces. 

1 2 



m. 



+ r(4 - ^cospf,) 



+ 0{1/N) 



{5i 



The smallest eigenvalue of this matrix is and corresponds to the p = matrix element. 
This in turn corresponds to the zero momentum modes of the quarks which from now on 
we will simply refer to as "zero modes" . 
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Figure 8. r = 0, mo = 0,A'^ = 2. The diamonds are the values of niq from the numerical 
simulation for Ix — S, It — 16. The sohd lines are the large N numbers with the zero mode 
included. From right to left they correspond to (L^, = 8,Lt — 16), (L^, = 16, = 16), 
{L^ = 32, Lt = 32), {L^ = 64, Lf = 64). The dotted lines are the large numbers 
with the zero mode excluded. From right to left they correspond to {L^ = 16, = 16), 
(La; = 8,Lt = 16). The solid vertical line denotes the infinite volume (3i^ from the large N 
calculation. For Pi < the model is in the broken phase. 




For small rUg the condition number c of this matrix is: 

c ~ — ^ for r = 



64r^ 

c — — for r — 1 (59) 



ml 
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A large condition number will make the inversion of M'^ M very slow. An important ob- 
servation can be made by noticing the dependence of the condition number on r. This 
suggests that performing the simulation with smaller r will yield a quite faster inversion. 
It is possible that this may also be the case for QCD. 



-10.0 




0.0 0.2 0.4 0.6 0.8 

Figure 9. r = 1, /9i = 2.5, N = 2, = 8, Lt = 16. The lines arc the large N effective 
potential with uiq = —0.2 (bottom) and mo = —0.5 (top). The crosses indicate the absolute 
minimum. Vgff is infinite at the spikes. 

But the unwelcomed effect of the zero modes on a finite lattice is not limited to large 
inversion times. Because on a finite lattice their effects are not suppressed by the measure 
but instead by an inverse volume factor, it turns out that in certain cases they severely 
obscure the physics. 
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Figure 10a. r — 0, mo = 0, N = 2, Lj. = 16, Lt — 16. The solid lines are the real 
part of the inverse a propagator for external 4-momentum q — {irricr, 0, 0, 0) at large 
with the zero modes included. From left to right they correspond to f3i = 2.3, 2.25, 2.2 or 
equivalently to o"s = 0.278, 0.323, 0.366. The functions have discontinuities that are denoted 
with the doted lines. The star denotes the relevant zero for the /3i = 2.2 line. 

In Figure 8 we plot niq = (a) vs. Pi for r — 0, uiq — 0, — 8, and = 16 (same 
as Figure 3). By simply looking at the numerical results we would not only be unable to 
estimate the critical point but also we would be unable to sec any indication of a phase 
transition. The large result on the same size lattice also has the same problems. As we 
increase the lattice size in the large N calculation (sohd lines from top to bottom) we see 
that a picture of an order parameter slowly materiahzes. At L^. — 64, = 64 a fairly good 
prediction of the large A^ infinite volume critical point is achieved. If we now do the same 
large A^ calculation but neglect from the momentum sums the zero modes, we obtain as a 
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result the two dotted lines for = 8, Lt = 16 and = 16, Lt = 16 (from left to right). 
We see that neglecting the zero modes on an = 8, = 16 lattice gives very similar 
results as the ones obtained on a L^. = 64, = 64 lattice with the zero modes included. 
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Figure 10b. r = 0, mo = 0, N = 2, jSi = 2.3 The lines are the real part of the inverse a 
propagator for external 4-momentum q = {im^, 0, 0, 0) at large N. The dotted line is for 
Lj, = 16, Lf = 16 with the zero mode included, the solid line for = 32, Lf = 32 with the 
zero mode included, the dot-dash line is for = 16, Lt = 16 with the zero mode excluded, 
and the dashed line is the infinite volume result from the asymptotic expansion. The stars 
denote the relevant zeros. 



If we plot the effective potential of eq. [T^ vs. ag for a finite lattice we will obtain 
a result as in Figure 9. The "spike" is a result of the presence of the zero mode on the 
lattice sum and extends to infinity. The presence of this spike may create thermalization 
problems if the initial configuration is chosen on the "wrong" side of the "spike" . Of course 
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the "width" of the spike is neghgible and therefore this problem may not be important. 

As mentioned at the end of the previous section, although one would expect to be able 
to measure the a mass in the r = case we were not able to do so. Large provides an 
explanation of this unexpected problem. In Figure 10a we plot the real part of the inverse 
o propagator, eq. 0, for a finite lattice with = 16, Lt = 16 and external momentum 
set to q = {irria, 0, 0, 0}. The sigma mass should be obtained at the zero of this function. 
We see that because of the presence of a discontinuity we do not obtain a root until m„ 
becomes heavy. The presence of this discontinuity is again due to the zero modes. 

In Figure 10b we plot the left most case of Figure 10a (dotted line). As we already 
mentioned there is no zero. If we increase the lattice size to = 32, Lt = 32 (solid line) we 
see that a zero develops. The infinite volume result from the asymptotic expansion (dashed 
line) has a zero nearby and no discontinuity since the zero modes are fully suppressed. If 
on the Lj, = 16, Lt = 16 lattice we now exclude the zero modes (dot-dash line) we see that 
the discontinuity disappears and a zero very close to the infinite volume result is obtained. 



4.3 The uItt = line on a finite lattice at large N 

The = zero line is of particular importance since it is there that the continuum chiral 
limit is obtained. As discussed in Section 3.2 this limit is obtained in the part of the m^r = 
line that corresponds to small quark mass. That region was presented in Figures 2a, 2b 
for —0.2 < TJiq < 0.3. If, for the same values of m^, the zero pion mass line is calculated 
on a finite volume one finds that the corresponding ranges of k and (3i change by almost 
an order of magnitude. To be more specific, for r = 1 the point of the the zero pion mass 
line corresponding, for example, to rUg = 0.015 has rriQ = —2.75 and f3i = 0.34 at infinite 
volume (see Figure 2a). The same point on an L^ = 8, Lt = 16 lattice has mo = —0.37 and 
/?! = 2.51. This change is fairly unusual and it may appear as if there is a contradiction 
between the finite and infinite volume results. In particular, in the infinite volume work of 
10| no point with m^ = was found for jSi > 1.41. To clarify this issue we plot in Figure 



11a the = line for an L^ = 8, Lt = 16 lattice for — oo < m„ < +oo. This line is 



calculated as follows: For a given value of m„ we use the finite volume version of eq. E6 



1 1 _|_ "'" w{k) ^ ^g^^ 



to calculate mo. Next, using this nig and mo we calculate Pi by combining equations [10 
and O The finite volume version of this is: 



2N 1 mg + rw{p) 
nig - mo L%Lt ^ g{p, nig) 

At point "1" of Figure 11a nig is very large and positive. As nig decreases, crosses zero 
and tends to very large negative values, we transverse the whole solid line from "1" to "2" 
to "3" ■ ■ ■ to "12" where nig is very large and negative. The "prongs" of this figure, points 
"3", "5", "7", "9", and "11", extend all the way to /?i = oo where they correspond to 
mg = 0, —2, —4, —6, and — 8, respectively. This singular behavior originates from the terms 
of the momentum sum in eq. ^ that correspond to the origins of the "Brillouin zones" 
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of the 16 species (see Table 1). When these terms are neglected from the calculation (the 
contribution of these terms disappears in the infinite volume limit) we obtain the dotted 
hne of Figure 11a. 




Figure 11a. r — 1^ N — 2^ — '^^ Lt — 16. The sohd line denotes the ttItt = line. 

To see more clearly how this behavior develops as we increase the volume we concentrate, 
as an example, on the "prong" that corresponds to ruq = since it is this "prong" that is 
most interesting in the recovery of continuum physics. A blown up picture of this "prong" 
is plotted in Figure lib. The outer solid lines correspond to the = 8, Lt = 16 lattice. 
Moving inwards the solid lines correspond to an = 16, Lt = 16 lattice and an = 32, 
Lt = 32 lattice. The dashed line corresponds to the L^ = 8, Lt = 16 lattice but with 
the singular terms removed. The dotted line is the infinite volume result of Figure 2a for 
—0.1 < niq < 0.1. Because the singularities are only quadratic they have disappeared in 
the infinite volume limit. The circle corresponds to the infinite volume nig — 0, 171^^ — 
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point. As the volume increases the whole "prong" is slowly mapped down to the infinite 
volume result of the dotted line. It is also very interesting to notice how close the infinite 
volume result is to the = 8, Lt = 16 lattice result with the singular terms removed. 
Similar behavior is expected for the other "prongs" as well. 
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Figure lib. Detail around the "prong" that corresponds to rUq = ("prong" "3" of figure 
11a). 



In 10 



it was shown that the = line is the phase line that separates the parity- 
fiavor symmetric phase ((|7r|) = 0) with the spontaneously broken parity-fiavor symmetry 
phase ((|7r|) ^ 0). We have confirmed that this is also the case on the finite lattice. The 
region inside the octopus-like graph of Figure 11a has ((|7r|) ^ 0). 

Finally, notice that the numerical results for m^r of Figure 7a indicate that on a finite 
lattice can be made quite small. The place where this happens is predicted by large 
quite well and it corresponds to a point on "prong" "3" that is quite far from where the 
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corresponding point of the infinite volume would be (somewhere on the dotted line on the 
"inside" of the "prong" ) . 

5 Conclusion and Summary 

In this paper we have found that the lattice version of the NJL model is an excellent toy 
model to investigate issues related to lattice QCD. We have used the large N approximation 
to leading order in to obtain non perturbative analytical results over almost the whole 
parameter range. By using numerical simulations of the model we estimated that for most 
of the quantities we considered the 1/N corrections are small. 
The main results of our investigation are listed below: 

1) With Wilson fermions we obtain at large N analytical expressions of the pion mass {ttIt,) 
and quark mass (mq) in lattice units as functions of the bare parameters of the model. We 
are then able to make exact statements regarding the approach to the continuum chiral 
limit. The "phase diagram" is presented in Figure 2. This may provide an insight on how 
the retrieval of the continuum chiral limit is achieved in QCD. 

2) At large and for Wilson fermions the a particle has mass proportional to the cutoff. 
Our analysis traces this fact to two related reasons. First, although the Wilson term 
has raised the masses of the doublers to the cutoff, it has not decoupled them from the 
theory. Through vacuum polarization these contribute to the a self energy and raise its 
mass. Second, although the Wilson term has not altered the low frequency behavior of 
the propagating quark, it has however altered its high frequency behavior. Again through 
vacuum polarization the high frequency modes contribute to the a self energy and also 
raise its mass. Such a phenomenon may also be responsible for the difficulty in observing 
a a particle in numerical simulations of QCD with dynamical Wilson fermions. 

3) The numerical simulation is performed on finite lattices. For naive fermions one would 
expect to be able to see some indication of the chiral phase transition as well as a cr particle. 
However, by simply looking at the graph of the vacuum expectation value vs. the coupling 
(see fig. 8) on an 8^ x 16 lattice one can not see any indication of a phase transition. Also 
the (7 particle in a 16^ lattice is either non existent or too heavy to be measured. Both of 
these unexpected results can be explained at large A^. The reason is traced to the existence 
of zero quark momentum modes that on a finite lattice are not sufficiently suppressed. The 
zero modes besides obscuring some of the physics are also probably partially responsible for 
the large inversion times in the HMC algorithm. To leading order at large A^ the smallest 
eigenvalue of the matrix that has to be inverted is and corresponds to the zero quark 
momentum mode. For small m^, the condition number of the matrix is for r = and 
QAr'^/mfj for r = 1. An important observation can be made by noticing the dependence of 
the condition number on r. This suggests that performing the simulation with smaller r 
will yield a quite faster inversion. It is possible that this may also be the case for QCD. 

4) The observables measured in the numerical simulation (chiral condensate, vacuum ex- 
pectation value, pion wave function renormalization constant, pion mass) have values that 
are in good agreement with leading order large A^. This provides a quantitative prediction 
for the size of the 1/N corrections. In agreement with the large A^ predictions discussed 
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in 2 and 3 above, the a mass was very heavy to give a good signal and was not measured. 
Also, measurements of the sigma width were not performed, but, as it will be discussed in 
5 below, we expect the corrections to the width to be large. 

5) For naive fermions we calculate at large N and with = 140 MeV the a mass {M„), the 
a width {To) and the quark mass (Mg)in physical units as functions of the cutoff. By setting 
Mq = 310 MeV we find = 726 MeV, = 135 MeV, and A = vr/a = 1150 MeV. 

is consistent with phenomenological expectations and A is consistent with the expectation 
that the cutoff should be close and below the mass of the lightest glueball (1550 MeV). The 
width however is underestimated. The reason is traced to the fact that to leading order in 
large N the width receives contributions only from the quark bubble and not from the pion 
bubble because the pion bubble is of order 1/A^. Because the phase space available for the 
cr to decay to two quarks is much smaller than the phase space to decay to two pions the 
pion loop contribution although of order 1/A^ is probably more important than the quark 
loop contribution. 

6) The above result can also be used to make an interesting observation. If the Higgs sector 
is the low energy effective field theory of a NJL model with exactly the same parameters 
as the low energy QCD except for = and F^^ = 246 GeV, then we find M^- = 
1915 GeV. This corresponds to ma = 2 where one would expect very large deviations from 
the low energy behavior of scattering cross sections. Although we have not calculated these 
deviations the value of the width serves as an indication of their size. In a way, departure 
from low energy behavior will be signaled by an increasing width of the a to two quark 
decay. At rricr ~ 2 the width is already fairly large. 

There are some interesting issues relevant to lattice work that have not been considered 
in this paper. It would be important to calculate the three and four point vertices and 
therefore be able to calculate scattering amplitudes and their departure from low energy 
behavior as well as the 1/A^ corrections to the width. It would also be interesting to study 
the NJL model at finite temperature and investigate the finite temperature transition in 
connection with the approach to the continuum chiral limit. Finally it would be important 
to include vector meson couplings (see, for example, fill]) and confirm that for the case 
of Wilson fermions the vector meson masses scale appropriately and do not become of the 
order cutoff as the a particle does. 
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Appendix A 



Below we give the lattice constants, defined in eq. ^ and that appeared in our analysis 
of Wilson fermions. 
r = 1 

ao,i = 0.0854 

ai,i = 0.2347 

ai,2 = 0.0260 

02,2 = 0.0597 

a2,3 = 0.0086 

= 0.0165 

ri = -0.0119 



Although they were never used, we have also calculated the constant defined in eq. ^ 
and also few more of the a„,m's. 

i?3 = -0.0016 
04,4 = 0.0050 
a4 3 = 0.0470 



"6,4 

0.1 



0.0388 (63) 



ao,i = 0.5373 

ai,i = 2.0721 

ai,2 = 1.3881 

a2,2 = 5.8935 

a2,3 = 6.1127 

03,3 = 25.633 

n = -0.3818 



The asymptotic expansions of the integrals Jo,i("^(?) and Ji,i(mg) defined in eq. [33 to 
order mj, are: 



JoAmq) = ao,i - 2m,rai,2 + m^[4r^a2,3 + n + Si + silogm^] + 0(mg) 
JiAmq) = ai,i - 2m<;ra2,2 + mg[4r^a3,3 - 01,2] + O(m^) (65) 



The leading order term of the asymptotic expansion of ^0,2(^^9) defined in eq. ^ is: 

JoA^q) = -'^1 - si- si logm^ + 0{mq) (66) 
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